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For non-zero cosmological constant A, we show that the barotropic FRW cosmologies as worked 
out in the comoving time lead in the radiation-dominated case to scale factors of identical form as 
for the Chiellini dissipative scale factors in conformal time obtained recently by us in Phys. Lett. 
A 379 (2015) 882-887. This is due to the Ermakov equation which is obtained in this case. For zero 
cosmological constant, several textbook solutions are provided as particular cases of A yf 0. 
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1. Parametric solutions of the barotropic FRW equations for non-zero cosmological constant 

Barotropic fluids continue to be of interest at the frontier of cosmology. In the past few years, they have been used 
as illustrative minimal models in the challenging problem of the dark sector of the Universe [l[, and currently they 
are considered as a viable partial solution for the late acceleration of the Universe due to their rapid fading to a de 
Sitter phase 0 . 

In this Letter, we are concerned with the solutions of the Einstein-Friedmann dynamical equations of barotropic 
FRW cosmologies with a cosmological constant A for the scale factor of the universe, a{t), with ' = d/dt 
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where k is the curvature of the universe, which is 0,±1, if the universe is flat, closed or open, respectively, p, p are 
the energy density and pressure, t is the comoving time, G is the universal gravitational constant, A is Einstein’s 
cosmological constant which is related to the pressure of the vacuum, and H is the Hubble expansion parameter. 
Since for a barotropic fluid the equation of state is 


P= ( 7 - 1)P, 

then the FRW equations can be decoupled and yield to 

ad -\- -b 7/t = — (7 + l)a^, 
o 

where 7 = I7 — 1 0. 

Let us introduce the transformation 


then z satisfies the linear equation 
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By one quadrature with positive initial conditions we have 
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z(a) = Cia — if j ^—1 

z(a) = €20^ — K, if 7 = —1 

and by using equation ([1]) we obtain the following parametric solutions 
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In addition, since z{a) is known from we substitute it back into and using ((3|) we obtain the equation of 
motion for the scale parameter in which the curvature of the universe does not occur explicitly 


d — C 2 a = 0, 

2. The radiation-dominated case, 7 = 1 


d-fa =-Ci 7 a-( 27 +i)^ if 7^-! 


if 7 = — 1 . 
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Interestingly enough, for the universe dominated by radiation one may recognize ® as the following Ermakov’s 
equation 


-A ^ _3 

a - a = —Gio . 
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From this equation, by multiplication with a and one integration, one can obtain 

a da 


t-to = 
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which is the solution given in the first line of ® in the case of radiation. 


K = ± 1 . The solutions of the above integral depend on the sign of the discriminant A = 
and are shown in Table HI 
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At = — 1 (top signs and pair of signs in the main square root) 


At = 1 (bottom signs and pair of signs in the 

main square root) 
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iy =pl ± VAcoshV 2 A(t — to) 
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TABLE I: Solutions for open and closed radiation-dominated universes, A 7^ 0 , A 



We notice that these solutions are similar to the recently introduced Chiellini dissipative scaling factors of 
barotropic FRW universes in conformal time because if we identify = A^, ci = —k, and 2 k = Ci in the 

Chiellini-damped solutions given in Q their discriminants turns into the discriminant in the non-zero cosmological 
constant cases discussed here. These analytical solutions are based on the Chiellini integrability condition for Abel 
^nations of the first kind [^, which also describe the (Wnamics of FRW cosmologies filled with bulk viscous fluids 
[3 [ 1 ] , as well as FRW universes filled with scalar fields Q ■ 


At = 0 . In this case, we have 
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re = 0 

7 = 1 


^sinh2A(t — to) :A>0 
^sin 2 X{t — to) : A < 0 


TABLE II: Solutions for the flat radiation-dominated universe, A 7 ^ 0, A 



The solutions of the above integral are given in Table HIl 

3. Cases of zero cosmological constant, A = 0 

Zero cosmological constant cases can be treated as particular cases of the formulas O- Let us denote the two 
integration constants from 0 as Cl = Oq^, C 2 = Oq^, which mean that all solutions pass through a(to) = ao, then 
the solutions of ([71) are 

a(t) =ao[l + (7 + l)^]^, if 7 7 Z -1 , k = 0 

t-tp _ 

a{t) = aoe “o , if 7 = —1, k = 0 

The first two solutions are obtained easily by integrating equations ([7|), and they represent the restricted form of 

1 

flat solution when to = for which a(t) = ao(^) j in the first case, while for to = aolnoo, we obtain the de 

t 

Sitter flat solution a(t) = e‘o in the second case. 

For the three textbook cases of 7 = 1 (radiation), 7 = ^ (dust), and 7 = — 1 (vacuum), by solving ([5]) with A = 0 
the solutions are presented in Tables m and IlYl 


7 = 1 

a{t) — ao^—1 + (1 + ^7^)^ 

7=1 

^A(t)-\/l -f A{t) — arccosh-\/l -I- A{t) — + %/2 — arccosh%/2 

7 = -1 

a{t) = ao sinh + arcsinhlj 


TABLE III: Solutions for k = -1 , A = 0, A{t) = 


7=1 

a(t) = ao^l-(l-f ^)2 

7=1 

AA{t)Al - A{t) + arcsin^l - A{t) = 

7 = -1 

a{t) = aocosh 


TABLE IV: Solutions for k = 1 , A = 0, A{t) = 

In the last case, 7 7 ^ 0, k 7 ^ 0, we will follow the substitution given by Lima see also Q 

u{a) = , (13) 

where A{t) = and differentiating u{a) with respect to t we have 

u = 2ju—. (14) 

a 

By inverting equation m which will give a and finding a from we obtain the first order equation that gives t(u) 
as the incomplete beta function 

^ - u)~^, 
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where /r = (7 + l )/27 and D = aoK^^^'*'/ 27 . We differentiate again with respect to u and obtain the equation 


, . drt 


(1 - m) - 
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dt 

du 


= 0 , 


(16) 


which is the hypergeometric equation with coefficients a = 0 , 5=i — /j, = —and c=l — /i = 5 — 
The general solution of (fTBl) is the linear combination 
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However, according to formula 15.4.1 in 0 , the first hypergeometric function is equal to unity. Thus, the general 
solution can be written in the form 


t{u)=a + Pu-^ + + 


1 1 
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It is convenient to pass to the complementary variable 1 — u in the above hypergeometric function using the formula 
15.3.6 in [ 1 ^. Then, one obtains 


t{u) = a + Py/^ Q + 
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The first hypergeometric function in (IT^ can be shown to be = which leads to the final result 

1 i-i_J- „ / 1 13 
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t{u) =to- 2/3(1 - u) = u 2 + 2’7 23^1 (1, ^ 1 - m 
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where to (the origin of time) depends on 7 

to = a + (dy/iT . ( 21 ) 

By choosing carefully the integration constants a and B in (I20L one can obtain the textbook solutions presented in 
Tables Uni and HVl 

The general solution (1201) is not valid when c is either equal to unity or any negative integer, i.e., when 7 = l/(2n+l), 
n = —1,1,2... . In these cases, the general solution has been given in Q, see also Q, and involves the logarithmic 
solutions of the hypergeometric equation as the other linearly independent solution. This is due to degeneracy for 
c = 1 and meaningless of one of the hypergeometric series in all the other negative integer cases of c. 
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